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I. INTRODUCTION 



In the view of relativistic quantum mechanics, a particle moving in a potential field is described 
particularly with the Klein-Gordon (KG) equation. Solutions of the one-dimensional KG equation 
have been received great attention for some potentials [1-2]. The relativistic quantum mechanical 
problems that can be solved exactly are very restricted [3]. In the present work, we obtain approx- 
imate analytical energy eigenvalues, normalized wave functions and scattering phase shifts for the 
Yukawa potential [4] 

V(r) = -^e~ ar . (1) 

where a is the screening parameter and 77 is the strength of the potential. The Yukawa potential has 
many applications in different areas of physics: in high-energy physics as a short-range potential 
[4], atomic and molecular physics as a screened Coulomb potential and plasma physics as the 
Debye-Hiickel potential [5]. In recent years, considerable efforts have also been made to study the 
approximate bound state solutions of the Yukawa potential in the non-relativistic domain by using 
different methods [5-9]. 

On the other hand, the position-dependent mass (PDM) formalism [10] has many applications 
in different areas, such as impurities in crystals [11], the study of quantum wells and quantum 
dots [12] and semiconductor heterostructures [13]. In recent years, the relativistic extension of 
the position-dependent mass formalism has been studied by many authors for different types of 
potentials [14-16]. 

The organization of this work is as follows. In Section II, we study the approximate bound 
state solutions and corresponding normalized wave functions for the Yukawa potential. We list 
some numerical results for the cases of PDM and constant mass presented in Table I and II. In 
Section III, we deal with the approximate scattering state solutions of the Yukawa potential and 
give analytical expressions for the phase shifts. In Section IV, we give our conclusions. 

II. BOUND STATE SOLUTIONS 

The radial part of the effective-mass KG equation is written as [17] 
d 2 <p{r) f £{£+!) 



- ^KWc 4 - (E 2 - 2EV(r) + V 2 (r))] J#r) = 0, (2) 



dr 2 \ r 2 h 2 t 

where I is the angular momentum quantum number, E is the energy of the particle and c is the 
velocity of the light. 



In recent years, the following approximation 

V^ W (i-e-^ r ( 3 ) 

has been used [18, 19] instead of the centrifugal term in the wave equations to obtain the solutions 
with any £ values. It has a good accuracy for small values of the potential parameter a [18, 19]. A 
remarkable approximation is proposed by Alhaidari [20] where the author suggested, for the first 
time, an approximation for the orbital term 1/r in the Dirac equation not for the 1/r 2 term. By 
using this approximation, it is possible to find the approximate bound state solutions of the Dirac 
equation for coupling to pure 1/r vector potentials with any I values for higher excitation levels 
with more accuracy than using the traditional approximation for the 1/r 2 term [20] . 
We define the mass function as 

m ( r ) = m ° + e 2ar _ 1 > ( 4 ) 

where mo and mi are two parameters and mo will correspond to the rest mass of the KG particle. 
Using the approximation given in Eq. (3) Yukawa potential becomes 

nr) = -2a VT ^^, (5) 

Inserting Eqs. (4) and (5) into Eq. (2) and taking a new variable z = (1 — e~ 2ar )~ 1 (z — > oo for 
r — > and z — > 1 for r — > oo), we obtain 

z{l-z)*P£ + (l-2z)^ 
dz z az 

/3 2 , 2 4 ^ 1 , f^momic 4 £/3 2 rA 1 



+ | w + 1) _|_ Wc< _ E2) _L_ + 

'/3 2 m 2 c 4 2 2 \ 1 - z 



2a 2 2a 



z 



where /3 2 = l/h 2 c 2 . Taking the form of the wave function 

( / ) (z) = z x ^l-z) x ^(z), (7) 
and inserting into Eq. (6), one gets a hypergeometric-type equation [21] 
z(l-z)^ + [l + 2A 1 -2(A 1 + A 2 + l)z] #(z) 



dz 2 ' dz 

2™2„4 



+ {-A 2 -A 2 -A 1 -A 2 -2AiA 2 + ^(^ + l) + ^i^-/3V}V(^)=0, (8) 



where 



^ = ^(^-^)-^ + ^ + ^-W, (9) 

/3 2 



A i = i^Kc 4 -^ 2 ). (10) 



Comparing Eq. (8) with the hypergeometric equation of the following form [21] 

z(i - z) y " + [& - (6 + 6 + 1 W - 66y = o , (ii) 

we find the solution of Eq. (8) as the hypergeometric function 

iKz)=2F 1 (Sufr,fr,z). (12) 

where 
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& = A 1 + A 2 + Ul + Ji + M{l + l) + ^4^-Af3W ), (13) 



6 = A: + A 2 + ± (l - ^1 + U{1 + 1) + - 4/3V ) , (14) 

& = l + 2Ai. (15) 

From Eq. (7), we obtain total wave function 

<f>(z) = Nz^(l - z) A2 2^(6,6,6,^) • (16) 

where N is normalization constant that will be derived in Appendix A. When either £i or £2 equals 
to a negative integer — n, the hypergeometric function ip(z) can be reduced to a finite solution. 
This gives us a polynomial of degree n in Eq. (12) and the following quantum condition 



Ai + A 2 + - + -\/l + 4^ + l) + ^-^-4/3V =~n (n = 0, 1, 2, . . .). (17) 



It is the relativistic energy eigenvalue equation for the Yukawa potential within the PDM for- 
malism. Defining two new parameters such as 



R2 m 2 4 

£(£) = y 1 + 4^(£ + 1) + - 40V , (18) 

1 = 4^2 K c ~ E ) + + A ^ ; A ^ = 2^ + ~ & V ' ( 9) 

we get the approximate energy eigenvalues as 

£(l) + 2n + 1 



x 



a 2 



[C{1) + 2n + 1]2 + ^2] _ + 2n + 1]2 - 4A(/3)] 2 I . (20) 



Table I presents the comparison of our numerical results for the case of constant mass (mi = 0) 
with the ones given in Ref. [22]. We restrict ourselves for only s-states and take mo = 1 because of 



the computation in Ref. [22]. Our parameters rj and a correspond to A and k(= coX), respectively. 
The relativistic energy is obtained as Er = rj exact (^s/l — A 2 — 1^ + 1 in Ref. [22]. The same 
numerical values of r] exact is used with Ref. [22] to compare our numerical results. We plot the 
1/r and the approximation 2ae~ ar (l — e~ ar ) versus r. It seems that the energy eigenvalues have 
a good accuracy up to the values of r\ < 0.25 and a < 0.30. Table II presents numerical energy 
eigenvalues for the case of position dependent mass including also for the case of constant mass 
with different values of (n,£). 

Now, let us study the results of our formalism for the case of the Coulomb potential. 



1. Relativistic- Coulomb Limit 



In the limiting case a — > 0, the Yukawa potential reduces to 

V{r) = -1, (21) 

which is an attractive Coulomb potential received great interest not only in the case of constant 
mass [23-25] but also within the position-dependent mass formalism [26]. 

In order to compare our results for the bound states with the ones obtained in Ref. [27], we 
expand the mass function, Eq. (4), into Taylor series 

, . mi mi , ^ Mi . . 

m(r) -> m --i + -^ + ... = M + ^ + ..., 22 

a->0 2 loir r 

The parameters b and tuq used in Ref. [27] are defined as b — > Mi , mo — > Mq . Taking the vector 
part of the potential which is equal to the scalar part as stated in Ref. [27] as V(r) = S(r) = — j 
and inserting the mass function into the KG equation including the scalar potential 

d 2 (j){r) \ 1(1 + 1) + 1 , 2 4 + 2m(r)5(r) + S 2^ _ ( E 2 _ 2EV{r) + V 2 {r))} \<P{r) = , 

he 2 J 

(23) 



dr 2 r 2 



gives the following equation 



dr 2 A t-A2--A 3 - 2 ^(r) = 0, (24) 
where 

Al = /3 2 (M^c 4 - E 2 ) ; A 2 = /3 2 (2M Mic 4 - 2M c 2 t] - 2Erj) ; A 3 = (3 2 (Mfc 4 - 2Mic 2 t] + f3£{£ + 1)) , 



(25) 
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The wave function is written as 



0( r ) = r K e- Air f{r) . 



(26) 



where we set k{k — 1) = A%. Thus, inserting Eq. (26) into Eq. (24) and using a new transformation 
z = 2A\r, we obtain 



^ + <*-^ + <-«-£)/w = .. 

which has the form of the Kummer differential equation [21] 

xy"(x) + (c - x)y'(x) - ay(x) = . 



So, the solution of Eq. (27) is given by 



A, 



In order to get a finite solution, the following condition must be satisfied 



2A-\ 



-n (n = 0,1,2,...) 



(27) 



(28) 



(29) 



(30) 



We get the bound state energy eigenvalues for the Coulomb potential as 

Con/. _ M C 2 



ft. 



4/3V + [iV + ^l + rf + 4f3 2 £(£+l) 



x <^ 4/3^(Mic^ - if) 



+ (N + y/l + rf + Af3 2 £(£ + 1) ) y |jV + Vl + + 4/3 2 ^ + 1) ] 2 - r/'(M lC 2 - 2t/) j . (31) 



where 



iV = 2ra + 1 and r/ = 4/3M lC 2 (M iC 2 - 2rj) 
This result is the same with the one for £ = given in Ref. [27]. 



(32) 



III. SCATTERING STATE SOLUTIONS 



Now we turn to the solution of the Eq. (2) to obtain the scattering states for the Yukawa 
potential. We use a new variable s = 1 — e~ 2ar (s — > for r — > and s — > 1 for r — > oo) and obtain 



+ 



1 ^(s) J<P(s) 
ds 2 ds 

^-l - — 2 (m omi c* - £(£ + 1)) - - - E>) — 



+ { - 



/3 2 mfc 4 \ 1 — s 
Act 2 " 



4>(s) = o. 



(33) 
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Defining the trial wave function 

4>{s) = s kl (l -s) k2 if>(s), (34) 
and substituting into Eq. (33), we obtain a hypergeometric-type equation for ip(s) 

s(l-s)^ + [2k 1 -(2k 1+ 2ik' 2 + l)s} d ^ s) 



ds 2 v * 1 ' 1 ds 

2_2 „4 T?o2„ 



+ { _ 2ihk' 2 -h- £(£ + 1) - + ^M«) = , (35) 



where 



4a 2 

The solution of Eq. (35) is a hypergeometric function 



where 



P = h + ik' 2 + \ - fPrfi -i-AE*- my) - + ^ , (39) 



'/? 2 m 2 c 4 
4a 2 


- /3 2 r/ 2 - 


4a 2 


(E 2 


— WqC 4 ) 


/3 2 mow,ic 4 
2a 2 


Ep 2 i] 
+ 2a 


'/? 2 m 2 c 4 
4a 2 


- /3 V - 


P 2 
4a 2 


(E 2 


— rrigC 4 ) 


2a 2 


E(3 2 r] 
2a 



or 



(36) 



fci = i { 1 + \J 1 + 4^ + 1) + ^-f^ - 4/32^2 | j 

k 2 = ik' 2 ; k' 2 = \R {E 2 -m 2 c*) . (37) 



Ms) = 2 Fi(p,q;r;s), (38) 



r = 2fei. (41) 
From Eqs. (34) and (38), we write the wave function of the scattering states 

cP( s ) = s ^(l- s y k 2 2 F 1 (p,q;r;s), (42) 



0( r ) = (1 - e -2orjfci e -2ifcSor g; r; j_ _ g -2ar^ _ (43) 

To obtain a finite solution, p or q must be a negative integer. This gives the following equality 



. //3 2 m 2 c 4 n0 (3 2 . . , /3 2 m mic 4 £/3 2 ?7 

(44) 

which is the same energy eigenvalue equation given in Eq. (17). We write the asymptotic form of 
the wave function given in Eq. (43) to obtain the scattering phase shifts. For this purpose, we use 
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the property of the hypergeometric functions [21] 

2 F 1 (a,b;c;y) = ^ C ^ C ~ ° ~ S 2^1(0, b; o + 6 - c + 1; 1 - y) 
1 (c — a)L (c — 0) 

+ r(C) r ( ( a a ) r ( fc ) ~ ^ ( X ~ ^) C ~ a " b 2^i(c - a, c - b; c - o - & + 1; 1 - y) , (45) 
and 2-fi(a, c; 0) = 1, we obtain the wave function for the limit of r — > 00 

</,(> _^ oo) ^ fl _ e -2arxfc! f r(2fci)r(-2zfc 2 ) -2ik' 2 ar 

n \ r(fei - *A4 - >t(fci, fe))r(fci - ^2 + ^1,^2)) 
, r(2fc 1 )r(2ifc 2 ) 2ik , ar \ 

T{k 1 +ik' 2 -A{k l ,k 2 ))T{k 1 +ik' 2 +A{k l ,k 2 )) J ' 1 j 



which could be written 



0(r^oo) - e- 2ar ) fci r(2A;i 



r(2ife 2 ) 



r(fci + ik' 2 - A(h, k 2 ))T(k 1 + ik' 2 + ^.(An, fc 2 )) _ 



^—2ik' 2 otr 



i 2ik' 2 ar l , , - , 

r(fci + ^2-^1,^2)^(^1 + ^2 + ^1,^2)) 



where 



■*(*i> ^) = V^dr- " W " £5 " -o^ 4 ) - 0I2 + -£r • ( 48 ) 



4a 2 4a 2 v 7 2a 2 2a 

From Eq. (47) we obtain 

T(2i£: 2 



<t>(r -> 00) -> 2(1 - e - 2ar ) fci r(2fei; 



r(fci + ife 2 - ^(fei, fe 2 ))r(fei + ik! 2 + ^(fci, fca)) 

sm ( 2a ^ r + 1 + ar5 r (fcl + rk< 2 - A(kJ 2 ml + ig + ^ ) ' (49) 

and, consequently, the phase shifts 5g are obtained as 

St = {i + l) l + J = ( ^ + l) l + flrg r (fcl + i% - ^ + ^k' 2 + ^ **)) • (50) 

IV. CONCLUSION 

We have studied the approximate bound and scattering state solutions of the effective mass 
KG equation for the Yukawa potential. We have obtained the energy eigenvalues, normalized wave 
functions and scattering phase shifts approximately as well as for the constant mass case. We have 
presented our numerical results of the energy eigenvalues in Tables I and II. We have compared 
them for the constant mass case with the ones obtained in the literature. We have also studied the 
bound state solutions of the Coulomb potential obtained from the limiting case of a — > with the 
position-dependent and constant mass cases. We have observed that the results obtained for the 
Coulomb potential are the same with the ones obtained in the literature. 
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Appendix A: Derivation of Normalization Constant 

By using Eq. (17), the wavefunction in Eq. (16) is written 

4>{z) = Nz Xl (1 - z) X2 2 F 1 (-n, n + 2Ai + 2A 2 + 1; 1 + 2Ai; z) , (Al) 

We use a new variable s = z~ l (s — > for z — > oo and s — > 1 for z — > 1) to normalize the 
wavefunction. For this purpose, we substitute s — > 



<£(iy^) = N'(l - s) Al (1 + s) A2 2 F 1 (-n, n + 2X 1+ 2A 2 + 1; 1 + 2A i; i^i) , (A2) 



where N' = N2-^ +x ^ (-l) 2Aa . 
Using the following equality [21] 

pW) {x) = ( 1 + Q )n 2Fl( _ n) n + a + /3 + 1;1 + Q . 1Z£) ) (A3 ) 



n! 



the hypergeometric functions in Eq. (A2) could be written in terms of the Jacobi polynomials 
pji*'^\x). Here, (/c) n is defined («) n = , and setting a = 2Ai and /? = 2A 2 in Eq. (A3), we 

rewrite Eq. (A2) 

H^) = N'(l - s)^(l + 8 ) Aa »! r( ^ ± f^ P i 2Al ' 2A2) (^) • (A4) 
where the Jacobi polynomials P^°'^(,z) are defined [21] 



pWM -if r(-n + fc)r(a + 6 + n + fc + l)r(a + n+l) 1 _ fc , 

r (" n ) r(a + 6 + n + l) T(a + A: + l)fc! 11 Xj ' (AbJ 

Using Eq. (A5) in Eq. (A4) we obtain the wavefunction 

0(i^) = AT'S(n, fc)(l - s) Al+fc (l + s) x * , (A6) 

where 

5V„w-rn+2A^V 1 r(-n + /c)r(2A 1 + 2A 2 + n + fc + i) l 

S(n, fc) - r(l + 2A0 ^ S ^7Z^ r(2A 1 + 2A 2 + n + l) r(2A x + k + 1) " (A7) 



9 



The normalization condition \4>(^2^)\ 2 ds = 1 gives 



|iV| 2 |S(n, k)\ 2 f\l - s) 2Xl+2k (l + s) 2X2 ds = 1 . (A8) 
Jo 



Comparing the last integral with the following [21] 

T(6)T(u - 5) 



Jo 



2 F 1 ( 7 ,S;u;z), (A9) 



r(iz) 

and setting S = 1, z = — 1, 7 = — 2A 2 and v = 2Ai + 2k + 2 we find the normalization constant as 

\Nf = a 2 ^ + 2fc + 1 . (A10 ) 

1 1 |S(n,/e)| 2 2 Fi(-2A 2 ,l;2Ai + 2A; + 2;-l) 
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TABLE I: Comparison of Klein-Gordon ground state energies for different parameter values. 





a 




Vexact 


Ref. [22] 




our results 


0.125 


0.01250 




0.83072460 





.993484 




0.998702 




0.06250 




0.30947218 


0. 


.997573 




0.999999 




0.09375 




0.12370738 





.999030 




0.999542 




0.12500 




0.02452195 


0. 


.999808 




0.998100 




0.14375 




0.00187260 


0. 


.999985 




0.996759 


0.25 


0.0250 




0.88881431 





.971776 




0.994130 




0.1250 




0.35655334 


0. 


.998678 




0.999960 




0.1875 




0.15650395 


0. 


.995030 




0.998556 




0.2500 




0.04068600 


0. 


.998708 




u.yyoioo 




0.3000 




0.00257940 


0. 


.999918 




0.985799 


TABLE II: Energy eigenvalues of the Yukawa potential for different values of n 


and £ (in 








h = mo = c = 


1 unit). 
























n t 


V 


a 


mi = 


mi = 0.1 




mi = 


mi = 0.1 





0.1 


0.01 


0.999181 


0.411464 




0.998173 


0.394898 




0.01 


0.1 


0.995475 


0.859773 




0.994475 


0.855513 


1 


0.1 


0.01 


0.999987 


0.614868 




0.998985 


0.602709 




0.01 


0.1 


0.980294 


0.900967 




0.979294 


0.898992 


1 


0.1 


0.01 


0.999911 


0.638787 




0.998910 


0.627267 




0.01 


0.1 


0.954438 


0.887484 




0.953438 


0.885983 


2 


0.1 


0.01 


0.999913 


0.712338 




0.998912 


0.702947 




0.01 


0.1 


0.954440 


0.884025 




0.953440 


0.882563 


1 


0.1 


0.01 


0.999622 


0.725851 




0.998622 


0.716879 




0.01 


0.1 


0.917015 


0.852051 




0.916015 


0.850766 


2 


0.1 


0.01 


0.999200 


0.748196 




0.998199 


0.739937 




0.01 


0.1 


0.866525 


0.801327 




0.865525 


0.800141 


10 


0.1 


0.01 


0.994432 


0.888765 




0.993432 


0.885794 


5 


0.1 


0.01 


0.987613 


0.897157 




0.986613 


0.894680 


10 


0.1 


0.01 


0.978199 


0.900957 




0.977199 


0.898987 
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FIG. 1: Comparison of 1/r (full line) with 2ae Qr (l — e ar ) for different values of a = 0.25, 0.30, 

respectively. 
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